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Abstract 

We study a diffusion model of phase field type, consisting of a system of two 
partial differential equations encoding the balances of microforces and microenergy; 
the two unknowns are the order parameter and the chemical potential. By a careful 
development of uniform estimates and the deduction of certain useful boundedness 
properties, we prove existence and uniqueness of a global-in-time smooth solution 
to the associated initial/boundary- value problem; moreover, we give a description 
of the relative w-limit set. 
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1 Problem setting 

The Cahn-Hilliard system: 

d t p- K Ap = 0, fi = -Ap + f(p), (1.1) 

describes diffusion-driven phase-segregation processes in a two-phase material body. Here 
p, with p(x,t) G [0,1], is an order parameter field interpreted as the scaled volumetric 
density of one of the two phases, k > is a mobility coefficient, and p is the chemical 
potential; f stands for the derivative of a double-well potential /. Customarily, the two 
equations (II. ip are combined so as to obtain the Cahn-Hilliard equation: 

d t p = KA{-Ap + f'(p)), (1.2) 

a nonlinear high-order parabolic PDE for the order parameter that has been studied 
estensively. With this procedure - we note for later reference - the chemical potential is 
left in the background; in particular, there is no need to take an a priori decision about 
its sign. 

To achieve their generalization of (II. 2p . Fried & Gurtin and Gurtin j8[ [10] propose: 
(i) to regard the second of (II. ip as a balance of microforces: 

div£ + 7r + 7 = 0, (1.3) 

where the distance microforce per unit volume is split into an internal part ir and an 
external part 7, and the contact microforce per unit area of a surface oriented by its 
normal n is measured by £ • n in terms of the microstress vector £0 (ii) to interpret the 
first equation as a balance law for the order parameter: 

d t p= -divh + a, (1.4) 

where the pair (h, a) is the inflow of p; (iii) to restrict the admissible constitutive choices 
for 7r,£, h, and the free energy density ip, to those consistent in the sense of Coleman & 
Noll [1] with an ad hoc version of the Second Law of continuum thermodynamics, namely 
a postulated "dissipation inequality that accomodates diffusion" : 

d t i) + (vr - p)d t p - £ ■ V{d t p) + • Vp < (1.5) 

(cf. eq. (3.6) in [10]). Within this framework, the following set of constitutive prescrip- 
tions is shown acceptable: 

V> = rj}{p, Vp), 

n{p,Vp,p) = p- d p ${p,Vp), (1.6) 
?(p,Vp) = «9 V p^(p,Vp), 

together with 

h = -MVp, with M = M(p, Vp, p, Vp); (1.7) 

: In [7], the microforce balance is stated under the form of a principle of virtual powers for microscopic 
motions. 
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moreover, it is shown that the tensor-valued mobility mapping M must satisfy the in- 
equality: 

Vp ■ M(p, Vp, p, Vp) Vp > 0. 
It follows from (TO]), (fOj) . fOD . and $TQ 1 that: 

<9 t p = div (MV (<9^(p, Vp) - div (dvpfto, Vp)) - 7)) + a 

(cf. eq. (3.17) in pH]); in particular, the Cahn-Hilliard equation (jl.2p is arrived at by 
taking: 

^(p,Vp) = /(p) + ^|Vp| 2 , M = «J, (1.8) 

and both the external distance microforce 7 and the order-parameter source term a iden- 
tically null. 

One of us proposed in [15] a modified version of Fried &: Gurtin's derivation, in which 
their step (i) is retained, but the order-parameter balance fll.4p and the dissipation in- 
equality (11 .51) are both dropped and replaced, respectively, by the microenergy balance 

d t e = e + w, e := — div h + a, w := —it d t p + £ ■ V \d t p) (1.9) 

and the microentropy imbalance 

dtf) > — div h + a, h := ph, a := per. (1-10) 

The salient new feature of this approach to phase-segregation modeling is that the mi- 
croentropy inflow (h,a) is deemed proportional to the microenergy inflow (h,a) through 
the chemical potential p, a positive field; consistently, the free energy is defined to be 

if; := e - p~\ (1.11) 

with chemical potential playing the same role as coldness in the deduction of the heat 
equation^] Combination of fll.9p - fll.lip gives: 

d t ip < -^(p" 1 ) + p" 1 h • Vp - 7T d tP + £ ■ V(d t p), (1.12) 

an inequality that replaces (jl.5p in restricting a la Coleman & Noll the possible consti- 
tutive choices. 

On taking all of the constitutive mappings delivering ir,£, p, and h, dependent in 
principle on p, Vp, p, Vp, and on choosing 

1> = $(p, Vp,p) = -pp + /(p) + ^|Vp| 2 , (1.13) 
compatibility with f ll,12p implies that we must have: 

n(p, Vp, = <9 P </>(p, Vp, p) = p - /'(p), 

?(p, Vp, p) = <9 Vp ^(p, Vp, p) = Vp, (1.14) 
^(P> Vp, = V 2 d^(p, Vp, p)= -p 2 p, 

2 As much as absolute temperature is a macroscopic measure of microscopic agitation, its inverse - the 
coldness - measures microscopic quiet; likewise, as argued in [15], chemical potential can be seen as a 
macroscopic measure of microscopic organization. 
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together with 

h(p, Vp, p, Vp) = —H(p, Vp, p, Vp)Vp, Vp ■ H(p, Vp, p, Vp) Vp > 0. 

If we now choose for H the simplest expression H = kI, implying a constant and 
isotropic mobility, and if we once again assume that the external distance microforce 7 
and the source a are null, then, with the use of (II. 14ft and (II. lip , the microforce balance 
(jl.3p and the energy balance (II .9p become, respectively, 

div(Vp)+p-f(p) = (1.15) 

and 

2pd t p + pd t p — nAp = 0, (1-16) 

a nonlinear system for the unknowns p and p that we supplement with homogeneous 
Neumann conditions at the body's boundary: 

d v p = d u p = (1.17) 

(here d u denotes the outward normal derivative), and with the initial conditions: 

p\t=o = Po, p\t=o = Po- (1-18) 

Needless to say, (11.151) is the same 'static' relation between p and p as (ll.ip 2 . Instead, 
(I1.16P is rather different from (ll.ip i . for a number of reasons: 

• (11.161) is nonlinear (whereas d t p — nAp = is a linear equation); 

• the time derivatives of p and p are both present in (jl.lfip ; 

• there are nonconstant factors in front of both dtp and dtp. 

Moreover, it should be possible to show that the initial/boundary- value problem (11.151) - 
(I1.18P has solutions p G [0, 1] and p > 0. 

We must confess that we boldly attacked this problem as is, prompted to optimism 
by the successful outcome of a previous joint research effort [5j |6], in which we tackled 
mathematically the system of Allen-Cahn type one arrives at via the approach in [15] for 
processes of phase segregation in the absence of diffusion. Unfortunately, system (11.151) - 
(I1.18P turned out to be too difficult for us. Therefore, we decided to study a regularized 
version of it, obtained by introducing two extra terms, e dtp in (I1.16P and 8 dtp in the 
left-hand side of (I1.15p . for small positive coefficients e and S. 

The introduction of the first term is motivated by the desire to have a strictly positive 
coefficient as a factor of dtp in (I1.16p . in order to guarantee the parabolic structure of this 
equation. As to the other term, on the one hand it gives (11.151) the form of an Allen-Cahn 
equation with right-hand side p; on the other hand, it assimilates our present model to 
the so-called viscous Cahn-Hilliard equations (see, e.g., [21 HH EES] and references therein). 
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With these measures, and taking k = 1 for simplicity, we write the following modified 
version of problem f ll.15p -f Tl.18p . with inversion of the order of the differential equations: 

{e + 2p)d t p + pd t p-Ap = infix(0,+oo), (1.19) 

5d t p - Ap + f (p) = p infix(0,+oo), (1.20) 

d„fj, = d v p = on T x (0, +oo), (1.21) 

p( • , 0) = p and p( • , 0) = p in f2, (1-22) 



where Q C M 3 is a bounded domain with a sufficiently smooth boundary T. We remark 
that such a regularized system has the typical features of a phase field model, but with 
a nonstandard equation (I1.19p for the chemical potential p, while quite often phase field 
systems feature temperature and order parameter as variables. 

By assuming, as we did in [HI IS], that /' is the sum of a strictly increasing C 1 function f[ 
with domain (0, 1) that is singular at the endpoints, and of a smooth bounded perturbation 
f 2 (to allow for a double- or multi-well potential /), we prove the existence of a strong 
solution (p,p) to (I1.19p - (ll.22p satisfying p > and < p < 1 almost everywhere in 
Q x (0, +00) (of course, the initial data have to meet the same requirements in Q). Our 
existence proof is rather standard; it is based on an approximation — > a priori estimates — > 
passage-to-the-limit procedure. Under additional assumptions, by using certain delicate 
iterative estimates, we also show that the component p is bounded above; this is probably 
the most difficult and technical part of the present paper. Boundedness of p is expedient 
to deduce that f'(p) is bounded as well; as a consequence, p stays away from the threshold 
values and 1. These boundedness properties are very useful in proving uniqueness of 
such solutions, since f'(p) can be treated as a Lipschitz-continuous function of p. 

As a final step, we deal with the long-time behavior of the system. We prove that each 
element (p^p^) of the u-limit set is a steady state solution of f ll.19p - fll.22p ; therefore, 
in particular, p u is a constant (cf. (I1.19p and fll.2ip ). This concludes our description of 
the contents of this paper. Needless to say, it would be interesting and challenging to 
study the singular limit of the solutions to fll.19p - fll.22p as e or 5 tends to zero, or both 
parameters do. We plan to undertake such a study in the near future. 

2 Main results 

In this section, we describe the mathematical problem under investigation, make our 
assumptions precise, and state our results. First of all, we assume O to be a bounded 
connected open set in IR 3 with smooth boundary T (to treat the lower- dimensional cases 
would only require minor changes). Moreover, for convenience we set: 

V := H l (Q), H := L 2 (fi), and W := {v G H 2 (n) : d u v = on T}, (2.1) 

and we endow these spaces with their standard norms, for which we use the self-explanato- 
ry notation || • ||y (but \\ -\\h denotes the norm of any power of H). We remark that the 
embeddings W C V C H are compact, because fl is bounded and smooth. Since V is 
dense in H, we can identify H with a subspace of V* in the usual way (i.e., so as to have 
that y*{u,v)y = (u,v)h for every u G H and v G V); the embedding H C V* is also 
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compact. As to the potential /, we assume that 



f = fi + f 2 , where functions f±, f 2 : (0, 1) — > R are such that 
fi is C 1 and convex, f 2 is C 2 , f 2 is bounded, 
lim f[(r) — — oo , and lim f[(r) — +oo. 

r\0 r/"l 



(2.2) 
(2.3) 
(2.4) 



For the initial data, we stipulate that 



Po £ V and p > a.e. in fi; 

PqEW, < p < 1 in fi; and /'(po) G 



(2.5) 
(2.6) 



We stress that the conditions in (I2.6P imply that 



(fi) and f(p )eH. 



(2.7) 



Indeed, W C C°(fi), assumptions ( 12. 3 j) hold, and, by convexity, — c < /i(po) < /i(l/2) + 
f[(Po){Po - I/ 2 ) for some c G R. 

Our aim is to solve problem fll . lQj) — fll .22j) in a strong sense, i.e., we want to find a pair 
(p, p) of such smooth functions satisfying suitable summability conditions and unilateral 
constraints that (jl.l9p - (jl.22p are made fully meaningful. Precisely, we fix a final time 
T > 0, we set Q := fi x (0, T), and we require that: 



Note that the boundary conditions (ll.2ip follow from (I2.8p - (I2.9I) . due to the definition 
of W in (12. ip . In conclusion, we look for (p,p) satisfying f!2.8j> — (I2.1ip and fulfilling the 

system 



Here is our main result. 

Theorem 2.1. Assume that (12.21) - 02 ,4p and (I2.5p - (l2.6p are satisfied. Then, there exists 
a pair (p, p) satisfing (12.81) - (12. lip and solving problem (12. 12p - (12.141) . 

Once existence is secured, one wonders about uniqueness. We are able to prove it for 
solutions having the following additional properties: 



peH\0,T- H)nL 2 (0,T-W), 

p E W 1,c °(0, T; H) n if^O, T; V) n L°°(0, T; W), 

p > a.e. in Q, 

< p < 1 a.e. in Q and f'(p) G L°°(0, T; 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



(e + 2p)d t p + pd t p - Ap = 

6d tP -Ap + f'(p)=p 
p(0) = p and p(0) = p 



a.e. in Q, 
a.e. in 
a.e. in fi. 



(2.12) 
(2.13) 
(2.14) 



p € L°°(Q); inf p > and supp < I. 



(2.15) 



Theorem 2.2. Assume that (j2"^?]) - (j2l| and (1231) - (I2TB1) are satisfied. Then, any two 
solutions to problem (12.121) - (12. 141) satisfing (|2TH]) - (I2.11I) and (12.151) coincide. 
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Interestingly, the additional boundedness conditions for (p, p) postulated above are ful- 
filled whenever the data of the problem have similar boundedness properties, in addition 
to (I23|)-(|23D. 

Theorem 2.3. Assume that (12. 2p - (12.61) and the following conditions are satisfied: 

Ho E L°°(£l); infp >0 and supp <l- ( 2 -16) 



Then, any pair (p,p) satisfing (12.81) - (12. lip and solving problem (I2.12l) - (l2.14p satisfies 
fl2~T5|) as well. 

Remark 2.4. Even though the regularity of the solution given by (I2.8l) - (l2.1ip and (12 . 1 5[) 
is completely satisfactory for our purposes, we observe that some further smoothness can 
be proved once the properties (12.151) are established. In that case, equation (I2.13P can 
be read in the form d t p — Ap = g with g G L°°(Q), whence further regularity for p can 
be derived, and a bootstrap procedure can start. Indeed, further regularity for p implies 
that stronger properties for p can be proved by (I2.12p . This improves the regularity of g 
and leads to an increase of the regularity of p. 

Once well-posedness on every finite time interval is ensured, one can study the long- 
time behavior of the solution. In particular, one can try to characterize the w-limit of any 
trajectory (/i, p) in some topology. We choose the weak topology of H x V and define 
such an w-limit as follows: 

weakly in H x V for some sequence t n / +oo}. (2-17) 

Our last result gives a relationship between such an cu-limit and the set of steady states, 
i.e., the set of the time-independent solutions (p s ,p s ) to (I2.12p - (l2.13p with homogeneous 
Neumann boundary condition satisfying natural regularity properties. Note that in such 
a case p s must be harmonic, thus constant, since f2 is connected. Therefore, a steady 
state is a pair (p s ,p s ) such that p s is a nonnegative constant and p s solves the following 
problem: 

p s G W, < p s < 1, f'(ps) G H, and - Ap s + f(p s ) = p s a.e. in Q (2.18) 
(there is no reason for p s to be constant, since / is not required to be convex). 

Theorem 2.5. Assume that conditions -d23D ; fl2"3D -fl27oD . and fl2~TB|) . are satisfied. 
Let (p, p) be the corresponding solution satisfing (I2.8p - (I2.11I) and (I2.15p . Then, the u-limit 
u(p,p) is nonempty, compact, and connected in the weak topology of H x V; moreover, 
each of its elements coincides with a steady state (p s , p s ) (that is to say, p s is a nonnegative 
constant and p s solves (12.181) ). 



Our paper is organized as follows. In the next section, we prove Theorem I2.1[ while 
Theorems 12.21 and 12.31 are proved in Section HI Our last section is devoted to the proof of 
Theorem 12.51 

Throughout the paper, we account for the well-known embedding V C L q (Q) for 
1 < q < 6 and the related Sobolev inequality: 

IMU<j(c7) < C\\v \\v for every v G V and 1 < q < 6, (2-19) 
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where C depends on Q only, since sharpness is not needed (the embedding V C L q (Q) 
is compact if q < 6). Furthermore, we repeatedly make use of the well-known Holder 
inequality, the interpolation inequality 

\\v\\ L r { a) < \\vf LP(n) IMfecn) for v e LP ^) n 

1 ■& 1-0 

where p, q, r <G 1, +oo , i? 6 0, 1 , and - = -H , (2.20) 

r p q 

and the elementary Young inequality 

ab < aa 2 + — 6 2 for every a, b > and a > 0. (2.21) 

Finally, throughout the paper we use a small-case italic c for different constants, that may 
only depend on Q, the final time T, the shape of /, the properties of the data involved in 
the statements at hand, and the coefficients e and 5; a notation like c a signals a constant 
that depends also on the parameter a. The reader should keep in mind that the meaning 
of c and c a might change from line to line and even in the same chain of inequalities, 
whereas those constants we need to refer to are always denoted by capital letters, just 
like C in (12~19D . 



3 Existence 

In this section, we prove Theorem 12. II Our method uses an approximation scheme based 
on a time delay in the right-hand side of (I2.13p . Namely, we define the translation 
operator T T : L 1 (0,T;iJ) — > L 1 (0,T; H) depending on a time step r > by setting, for 
v E L\Q,T;H) and for a.a. t E (0,T), 

(7 T v)(t) :=v(t-r) if t>r and (7 T v)(t) := fi if t < r, (3.1) 

and consider the problem obtained by replacing the right-hand side of (12.131) by T r /i, i.e., 
we look for a pair (/i r , p T ) such that 

(/i T , p T ) satisfies §ZM)-$rn$ (3.2) 

(e + 2p T )d t fi T — Ap T + p T d t p T = a.e. in Q (3.3) 

5d t p T — Ap T + f'(pr) = 7 T p T a.e. in Q (3.4) 

^ T (0) = po and p r (0) = po a.e. in Q. (3.5) 

For convenience, we allow r to take just discrete values, namely, r = T/N, where N is 
any positive integer. Our existence proof consists in two parts. Firstly, we check that 
problem (I3.2l) -( l3~5l) is well-posed (see the next lemma). Secondly, we let r tend to 0. This 
is done by proving a number of a priori estimates and using compactness and monotonicity 
arguments. 

Lemma 3.1. There exists a unique pair (p T ,p T ) solving problem (13.21) -f l3~5l . 

Proof. Recall that r = T/N. Hence, if we set t n := nr for n — 0, . . . , N, we see that 
problem (I3.2p -( l3"3|) becomes equivalent to a finite sequence of N problems that can be 
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solved step by step. However, instead of considering the natural time intervals [t n -i,t n }, 
n = 1, . . . , N, and glueing the solutions together, we solve N problems on the time 
intervals I n = [0, t n ], n = 1, . . ., N, by constructing the solution directly on the whole 



of I n at each step. These problems are the following: 

(e + 2p n )d t fi n - Ap n + {d t p n )Pn = and p n > a.e. in ft x (3.6) 

d u ^n(t)\v = for a.a. t G I n and p n {0) = p (3.7) 

< p n < 1 and 5d t p n - Ap„ + f'(p n ) = 7 T p n -i a.e. in ft x I n (3.8) 

dvPnit)\v = for a.a. t E I n and p n .(0) = po- (3.9) 



Their solutions are required to satisfy the regularity properties induction obtained by 
taking t n in place of T in (I2.8p - fl2.lip . The operator T r that appears on the right-hand 
side of (13. 8p acts on functions that are not defined in the whole of (0,T). However, its 
meaning is still given by (13.11) if n > 1, while we simply set 7 T p n -\ = p,Q if n — 1. 

Clearly, the solution (/i T , p T ) we are looking for is simply given by (/in, pn)- The above 
problems can be solved inductively, because the right-hand side of (13.81) is known at each 
step in the next lemma: one first solves problem fl3.8l) - fl3.9p for p n , and then problem 
fl3.6l) - fl3.7p for p, n . We note that the former problem is quite standard; the latter is a 
regular linear parabolic problem (the coefficient of dtp. n is > e) provided that dtp n is 
sufficiently smooth. That the inequality pi n > holds is not obvious. The uniqueness of a 
solution {pi n , p n ) satisfying smoothness properties analogous to fl2.8p - fl2.lll) is clear, and 
the existence of a variational solution is expected. However, not even the desired regularity 
is obviously guaranteed. Therefore, we provide a few arguments in this direction. 

Proceeding at an as-low-as-possible level of formality, we introduce a problem depend- 
ing on a positive parameter A and approximating problem fl3.8l) - fl3.9p . To begin with, we 
regularize j\ and j'2 by constructing certain suitable C 2 approximations /i,a,/2.a having 
bounded first and second derivatives. Precisely, we assume that f% x is bounded uniformly 
with respect to A; moreover, on thinking of f[ as a maximal monotone graph in 1 x 1, we 
assume that fi \ is convex and that f[ A is similar to the Yosida regularization of /{ (see, 
e.g., 0, p. 28]), in order to preserve the main properties of the latter (such a regularization 
is detailed, e.g., in [91 Section 3]). Finally, we set fx = fi t \ + /^.a- The approximating 
problem is: 

5d tPn - Ap n + f x {p x n ) = 7 Tf i n -i a.e. in ft x I n , (3.10) 
du(£(t)\ r = for a.a. t e I n , and p£(0) = Po ; (3.11) 

it has a unique smooth solution, which satisfies sufficiently strong a priori estimates to 
allow letting A tend to zero in (I3.10p - fl3.lip . This leads to a solution p n to problem (13.81) - 
(13.91) . which can be used to solve problem fl3.6p - fl3.7l) . Needless to say, the desired regu- 
larity for p n will follow once we prove suitable estimates uniformly with respect to A. We 
confine ourselves to derive the highest-order estimate, the others being quite standard. 

With a view toward assembling a proof by induction, we assume that 

^e^tU^jnrfUl') and /i„_i>0 for n > 1 (3.12) 

and we prove that 

\\pn\\w^(I n ;H)nm(I n ;V)nL^(I n ;W) + || f[,\(Pn) 1 1 < C r , (3.13) 

PneH 1 (In,H)nL™(I n ,V)nL 2 (I n ,W) and p n > (3.14) 
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(as anticipated in closing Section 2, in ( 13 . 1 3|) as well as in the following the symbol c r 
stands for one or another of a list of different constants that do not depend on A, but 
are allowed to depend on r). We remark that the induction procedure can actually start, 
because 7 T p n _i = p if n = 1 and, moreover, properties (12.51) and (12. 6p for p and p are 
fulfilled; these properties are also used at each step. 

We omit stressing the dependences on n and A, and write simply u and u Q for, re- 
spectively, d t p* and d t p^(0). By differentiating (13. lOj) with respect to time, we see that u 
solves the equation: 

8d t u -Au + fl x (p x n ) u = 9 t (T r // n _!) - fl x (p x n ) u a.e. in ft x I m (3.15) 

and satisfies both the Cauchy condition u(0) = uq and homogeneous Neumann boundary 
condition. Hence, by testing (13.151) by u and using the convexity of fi t \, we immediately 
obtain for t G I n that 

U<t)\\ 2 H+f 7 |V M | 2 <^|MWl+SUp|^ A |) ff u'+Wd^Pn^Wl^y (3.16) 

z Jo Jn z Jo Jn 

Now, we observe that the last norm is finite, in view of our assumption (13 . 1 2 j) . and that 
I/^'aI — c - Moreover, due to (I3.10p . we have that Suo = po + Ap — f'\(po)- Hence, uq is 
bounded in H, by (12 . 5 f) — ( |276T) and our choice of the approximation f\ of /. Therefore, 
thanks to the Gronwall lemma, we obtain: 

\\u\\L°°(i n ;H)nL 2 (i n ;V) < c T , whence ||p n ||wi I o°(7 n; ij-)nH 1 (/ n ;V) < c r . 

Next, coming back to (I3.10p . we deduce that —Ap x + f[ x {Pn) * s bounded in L°°(I n ;H) 
and hence, by a standard argument (for instance, by testing (I3.10P by / 1A (p^)), that each 
of —Ap x and f[xi.Pn) * s bounded. With this, given that the IV-estimate follows from 
elliptic theory, (13 . 131) is established, and we can let A tend to zero. We obtain: 

p n eW 1 '°°(I n ,H)nH 1 (I n ,V)nL c °(I n ,W), < p n < 1, and f[(p n )eL°°(I n ;H). 

At this point, we should prove (13. 14j) . However, we confine ourselves to derive a formal es- 
timate that clearly shows that the desired regularity for be deduced by regularizing 
the linear problem (I3.6p - (l3.7p (if the coefficient d t p n is replaced by a smooth function and 
the initial datum is regularized, by the same token d t p n is an admissible test function). 
For convenience, we write (13.61) in the form: 

(e + 2p n )d t p n + p n - Ap n = (1 - d t p n )p n ; 

next, we multiply this relation by d t p n and use the result in the calculation given below. 
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Since p n > 0, we find, for t G /„,, 



t r i rt r i 

2 i 1 II,. /jAI|2 ^ I I , r. \ I Q 1 2 i ll,, / + \||2 



/ l^r+olW*)llv < / / (e + 2pn)|fy*«r + dM*)llv 
jo z Jo Jn A 

= [[ {e + 2 Pn )\d t p n \ 2 + \\\p \\ 2 v + \ / / d t (\p n \ 2 + \Vp n \ 2 ) 
Jo Jn i z J j n 

= 7; WVoWv + / / ((e + 2p n )\d t p n \ 2 + p n d t p n + Vp n -Vd t p n ) 
1 Jo Jn 

Vo|| 2 /+ / / (1- d t p n ) Undt/In 



2 



o Jn 

t 



o 



< ^IIa*o|Iv-+ / II 1 + \dtPn(s)\\\L*(n) ||A*n(*)IU*(n) ||^/i n (s)IU2 ( n) 

1 ,, ,12 

< 2 WtoWv 



~ / / |^/i n | 2 + ^- / ||1 + |%>„(s)|||^ ||/x n (s)||^ds, 
^ Jo JO 25 Jo 



by the Sobolev and Young inequalities (12. 19ft and (I2.2ip . Then, the Gronwall lemma 
yields that 

\\dtPn\\L 2 (I n ;H) + \\ Pn \\ L°°(I n ;V) < Cm, (3-17) 

where M is a constant satisfying M > ||«o||v + \\dtp n \\L 2 (i n ;V)- By comparison in (13.81) . 
even Ap n is estimated in L 2 (I n ;H), since a bound for p n d t p n in the same space follows 
from (I3.17p . By elliptic regularity, we derive the desired estimate for p n in L 2 (I n ; W). So, 
the first assertion in (13.141) is established; it remains for us to show that p n > 0. This is 
done by testing (13.61) by —p~. We obtain, for t G I n , that 



\ f! d t ((e + 2 Pn )\p-\ 2 )+ fl \Vp-\ 2 
z jo Jn Jo Jn 

((£ + 2p n )d t p n {-p~) + {d t p n )p n {-p~) + Vp n ■ V(-/i„)~) = 0. 



'o Jn 

As p n > and p > 0, we deduce that 



/i~W| 2 < / {e + 2 Pn {t))\p n -{t)\ 2 < / (£ + 2p )|/io| 2 = 0, 
Jn Jn 

whence it immediately follows that p~ = 0, i.e., that p n > 0. Thus, the lemma is proved. 
□ 

Now that the well-posedness of problem (I3.2l) - (l3.5p is established, we perform a num- 
ber of a priori estimates of its solution. These estimates allow us to let r tend to zero, so 
as to prove our existence result for problem (I2.12l) -( [2.14l) . In order to make the formulas 
to come more readable, we shall omit the index r in the calculations, waiting for writing 
(p T ,p T ) only when each estimate is established. 

First a priori estimate. We observe that d t {{e/2)p 2 + p/i 2 ) ={{e + 2p)d t p + pd t p)p. 
Thus, testing (13.31) by p and integrating, we obtain, for t G (0,T), that 

/ (^ 2 + PV 2 )(t)+ ft |Vp| 2 = f (^pl + popl) =c (3.18) 



12 



Well-posedness and long-time for a nonstandard diffusion model 



This implies that 

||^r|U°°(0,T;iJ)nL 2 (0,T;V0 < c - (3.19) 

Second a priori estimate. This standard estimate for phase field equations can be 
derived by testing (13.41) by dtp. We get: 

l|Pr|| J f/ 1 (0,T;H)nL oo (0,T;y) + || / (Pt ) || L°°(0,T;L 1 (fi)) < C. (3.20) 

Third a priori estimate. We rewrite (13.41) as 

- Ap + f[(p) = -Sd t p - f 2 {p) + 0>, (3.21) 

and notice that the right-hand side is bounded in L 2 (0,T; H). Then, by applying a 
standard procedure (for instance, testing by f[(p)), and counting on elliptic regularity, 
we deduce that 

\\Pt\\l*{Q,T;W) + \\fl{pr)\W{0,T;H) < C. (3.22) 

Fourth a priori estimate. To derive the next inequality, we prefer to proceed formally, 
avoiding the A-regularization we used in the proof of Lemma [37X1 As for (13.161) . we obtain 
the following estimate: 

\\\dtp{t)\\H + fj^d tP \ 2 

<^||Ap -/((po)+/io||| + sup|^ A | ft \d t p\ 2 + [ [ (d t 7 T p) dtp. (3.23) 
z Jo Jn Jo Jn 

Once this inequality is established, our procedure is rigorous. The estimate of the last 
term requires now more care than before, because we aim to obtain bounds that are 
uniform with respect to r. We have: 

n(d t 7 T p) d t p= / d t pi(s - t) d t p(s) ds = / d t p(s) d t p(s + r) ds, 

a Jt Jn Jo Jn 

and we compute d t p from (13. 6p . On recalling that p > 0, we can continue as follows: 
J J (d t 7 T p) d t p = J J ( - ^ (A/i - p dtp)^j(s) d t p{s + r) ds 

- d t p(s) p(s) d t p(s + t) } ds 

e + 2p{s) J 

\^d tP \ 2 + 4p\\h {0>T]V) 

4 Jo Jn 

rt—T 

+ c I \\d tP (s + r)|| L 4 (n) ||V^(a)||H \\Vp(s)\\ LHn) ds 

\\d t p(s + r)\\ LHn) \\p(s)\\ LHn) \\d t p(s)\\ H ds. (3.24) 



rt—r 

+ c / 
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We need to estimate the last two integrals. As to the first, we begin by using the Sobolev 
inequality (I2.19p and the elementary Young inequality 02.21 j) . We find that 



t-T 



\\d tP (s + r)\\ LHn) \\Vfi(s)\\ H \\Vp(s)\\ L 4 {a) ds 
<\f T \\d t p(s + r)f v ds + cj)\Vp,(s)\\l\\Vp(s)\\lds 

+ c fMsWviWpisWv+WApis)^) ds, 



the last inequality holding because, thanks to elliptic regularity, \\v\\w < c(\\v\\v + \\Av\\h) 
for any v G V such that Aw G H and d v v\r — 0. At this point, we recall that p is bounded 
in H 1 (0,T;H)nL co (0,T;V), and [i in L 2 (0,T;V), by fEOUD and (I3TT9D . Moreover, we 
notice that 03.211) entails (formally, by testing it by — Ap(s)): 

\\Ap(s)\\l < 5^d t p(s)\\ 2 H + c(l + \\7Ms)\\ 2 H ) 

< ^\\dtp(s)\\ 2 H + c for a. a. s G (0, T), 

the last inequality being a consequence of (13. 19ft . Therefore, we can infer that 

t-T 

\\d t p(s + r)\\ LHn) \\Vp(s)\\ H \\Vp(s)\\ LHn) ds 
I fj a |V^p| 2 + c + c]\pl{s)\\1 \\d t p{s)f H ds. (3.25) 



< 



t-T 



'o Jn 

ft 



Passing now to estimate the last integral in (I3.24p . we have that 
\\9 t p{s + t)\\ L 4 (q) \\p(s)\\ L 4 (n) \\d t p(s)\\ H ds 

<~ f T ii^+^)ii 2 v^+^V(«)irvii^)iiH^ 

|V^p| 2 + \ \\d t p\\l^, T . tH) +cj)Hs)\\ 2 v \\dtP(s)\\ 2 H ds 

I H \Vd t p\ 2 + c + c tlHsWyWdtpisnlds. 
Jo Jn Jo 

With this and fl3T25|) . we see that fISTMj) yields: 

ft (d t 7 Tf ,)d tP <c+l f! |v<v| 2 + c f\\pt{s)\\l\\d tP {s)\\lds, 

Jo Jn * Jo Jn Jo 

so that ( 13.231) takes the form: 

^\\dtp(t)\\ 2 H + \ n \Vd tP \ 2 <c + c D^Wl \\d t p(s)\\lds. 



1 

< - 



< 



o Jn 
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Well-posedness and long-time for a nonstandard diffusion model 



Since p has been estimated in L 2 (0,T; V), the Gronwall lemma can be applied, so as to 
have that 

\\dtp T \\L°°(0,T;H)nL 2 (0,T;V) < c - (3.26) 

Finally, the same argument as in the derivation of ( 13 .22 j) yields that 

||Pt||l°°(0,7W) + ||/l(Pr)||L°°(0,T;H) < C. (3.27) 

Fifth a priori estimate. We formally test ( I3.3P by d t p, and obtain: 

e ff \d t p\ 2 + \ j \Vp{t)\ 2 
Jo J n z Jn 

< i II V/xo||lr - / / <V/i^ 

~ C+ h J ll^ 5 )l^ 4 ( n ) H/ i ( s )lli 4 (n) rfs + \J \\ d tP{s)\\ 2 H ds 

< C + C ^||^p s ( S )|| 2 y (||M S )||^+||V/i(^lli) eZ S + |^W(*)lll^ 

Jo 1 Jo 

where the last inequality follows from (13.261) and (I3.19p . Using (I3.26P once more, we can 
apply the Gronwall lemma and conclude that 

||AM|j* 1 (0,T;H)nZ, oo (0,T;V) < c. (3.28) 

Sixth a priori estimate. Recalling that < p < 1, and using the Sobolev inequal- 
ity fl2H, we get: 

|| (e + 2p)d t p + p d t p\\ L 2 [0yT . H) 

< (e + 2)\\d t fl\\L2(p,T;H) + ||p||l°°(0,T;L 4 (Q)) || <9tP|| L 2 (0,T;L 4 (n)) 

< c{\\dtP\\L2(0,T;H) + \\p\\l^{0,T;V) \\d t p\\l?{0,T;V)) ■ 

Since the right-hand side is bounded by (13.281) and (13.261) . a comparison in fl3.3[) shows 
that Ap is bounded in L 2 (0,T; H). Consequently, by elliptic regularity we deduce that 

\\Pt\\l2(0,T;W) < c (3.29) 

Conclusion. Collecting all the estimates we have proved, we see that 

p T — V p weakly star in H l (Q, T\ H) D L°°(0, T; V) PI L 2 (0, T; W), 
p T -»■ p weakly star in jy 1,oo (0, T; iJ) fl -ff^O, T; U) D L°°(0, T; W), 
f'i(pr) -> £ weakly star in L°°(0, T; if), 

at least for some subsequence \ 0. Thanks to the Aubin-Lions lemma (cf. [131 Thm. 5.1, 
p. 58]) and to similar results to be found in [13 Sect. 8, Cor. 4], we also deduce the 
following strong convergences: 

p T ^p strongly in C°([0, T]; H) H L 2 (0, T; 1/) 
p r p strongly in C° ( [0 , T] ; V) . 
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In particular, having recourse to a well-known monotonicity technique (see, e.g., [TJ 
Lemma 1.3, p. 42]), we conclude that < p < 1 and £ = f[(p) a.e. in Q. The 
strong convergence shown above also entails that /^(p-r) converges to /^(p), e.g., strongly 
in C°([0, T];H) (because f 2 is Lipschitz continuous), and that 7 T p T converges to p, e.g., 
strongly in L 2 (0, T; H). Finally, a combination of the above weak and strong convergence 
results with the Holder and Sobolev inequalities yields that 

p T d t p T — » p d t p weakly in /^(O, T; H), 
p T d t p T pd t fu weakly in L 2 (0, T; L 3 / 2 (fi)). 

Indeed, p T — ?■ p strongly in L 2 (0,T; L 4 (f2)), <9ip r — >■ <9 t p weakly in L 2 (0,T; L 4 (f2)), p r — p 
strongly in C*°([0, T]; L 6 (fi)), and d t p T <9 t p weakly in L 2 (0, T; L 2 (fi)). Therefore, it 
is straightforward to conclude that the pair (p, p) is a solution to problem fl2.12p -f l2.14p 
having the desired regularity f l2.8p -f l2TTTj) . that is to say, Theorem 12.11 is proved. 



4 Uniqueness and boundedness 

In this section, we prove Theorem 12.21 and Theorem 12.31 We first show our uniqueness 
result. 

Proof of Theorem 12.21 We take two solutions to problem fl2.12p - fl2.14p satisfying 
(12.1 5p in addition to (I2.8p - fl2.lip and label their components with the subscripts 1 and 2; 
in the following, the values of constants c may depend on these solutions. Moreover, we 
choose constants such M > and r*, r* G (0, 1) that p < M a.e. in Q and r* < pi < r* 
a.e. in Q, for i = 1,2. Finally, we denote by L the Lipschitz constant of the function 
t i — y t — fir), r 6 [r*,r*]. Having done this, we write (I2.12p for both solutions and take 
the difference. Then, we set p := pi — p 2 and p := pi — p 2 , test the resulting equality 
by p, and integrate, using the boundary condition. Due to the identity: 

{(diP^pi + 1p x d t p\ - (d t p 2 )p2 - 2p 2 dtp 2 } p = d t (pip 2 ) + 2(d t p 2 )pp + p 2 (d t p) p, (4.1) 
we obtain: 

/ ( £ 7; + Pi(t))\p(t)\'+ ft |Vp| 2 = -2 ff (d t p 2 )pp- ft p 2 (d tP )p. (4.2) 
Jn XA ' JoJn JoJn JoJn 

Furthermore, we write (12.131) for both solutions, take the difference, test the resulting 
equality by d t p, and add pd t p to both sides, for convenience. Then, we integrate, using 
the boundary condition, and easily obtain that 

^ ff \dtp\ 2 + h\ P (t)\\v= ff ((Pi-f(Pi))-(P2-f(p 2 ))+p)<V 
Jo Jn z Jo Jn 

<L f f \p\\d tP \+ f f \p\\d tP \. (4.3) 
Jo Jn Jo Jn 
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Now, adding ( 14. 2 p and (14. 3p and taking into account that pi is nonnegative, we get 



\ [ Ht)\ 2 + [ [\Vp\ 2 + 5 ft \d tP \ 2 + \\\p{t)\\v 
z Jn Jo Jn Jo Jn z 

<2 \\d t p 2 (s)\\ H ||p(s)|| £ 4(n) \\p(s)\\ LH n)ds+ / ||// 2 (s)|U»(n) \\dtp(s)\\ H \\p(s)\\ H ds 



+ L [ \\p(s)\\ H \\d t p(s)\\ H ds+ f \\p(s)\\ H \\d t p(s)\\ H ds. (4.4) 
Jo Jo 



To estimate the first integral on the right-hand side, we use the Sobolev inequality (1 2 . 1 9 [) 
with q = 4 and C the Sobolev constant, and we invoke the elementary Young inequal- 
ity (121TLD to obtain that 



\\dMs)\\ H \\p(s)\\ LH n) \Hs)\\L H n) ds < CM ||^ 2 ( S )|| H \\p(s)\\ v \\p(s)\\ v ds 

Jo 

< \ j\\p{s)\\lds + ^ j\d m {s)f H \\p{s)\\lds 

= \ ff |V/i| 2 + c ff \p\ 2 + c f\\dMs)\\ 2 H \\ P (s)\\ 2 v ds. 
A Jo Jn Jo Jn Jo 

The remainder of the right-hand side of (I4.4p is estimated as follows: 

t 

IMs)IU-(fi) \\d t p{s)\\ H \\p{s)\\ H ds 

+ L f \\p(s)\\ H \\d t p(s)\\ H ds + f\\p(s)\\ H \\d t p(s)\\ H ds 
Jo Jo 

<^ ff M + c f(\Hs)\\ 2 H +\\p(s)\\ 2 v )ds. 
A Jo Jn Jo 

Combining these estimates with (14. 4p . we immediately get 



n 

t 



W)? + \ loin ' V/i ' 2 + t loin M + \ 



<c (l + \\d t p 2 (s)\\ 2 H )\\p(s)\\ 2 v ds + c \\p(s)\\ 2 H ds. 



Since the function s i— > \\d t p2(s)Wjj belongs to L 1 (0, T), we can apply the Gronwall lemma 
and deduce that both p and p vanish. Hence, the two solutions coincide. □ 

We now turn to proving our boundedness result. 

Proof of Theorem EH. Let (p, p) be any solution to problem f[278|) ([27TTT) and ([2~T2l) - 
(I2.14p whose initial data have, in addition to (I2.5l) - (l2.6p , the further properties (12.161) . 
We show that the boundedness claims in (I2.15P actually hold true. 

With a view to proving that p satisfies the specified lower bound, we set: 

Po '■= ll/^o|U°°(fi) = supess/x (x) ; (4.5) 

x&n 
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we take any real constant k such that k > p 1 ^; and we introduce the auxiliary function 
X k e L°°(Q) defined for t) G Q by the formula: 

Xk(x,t) = 1 if p(x,t) > fc, and t) = otherwise. 

Then, we test (12.121) by (/i — k) + and integrate over Vt x (0,t) for any t e (0,T). The 
result is: 



J Q J «/ 



12 

^p|(/i-£;)+| 2 - /7 <9 t p/i (//-&)+=- / / fc <9 iP (// - k) + . 
oJn JoJn JoJn 

Given that p is nonnegative, this equality and the Holder inequality with ad hoc exponents 
lead to: 

£ ^\Mt)-k) + \\ 2 H + ff \V{p-k) + \ 2 
* Jo Jn 

<k \\Xk(s)\\ L 7/2 [Q) \\d t p(s)\\ L i4/ Hn) ||(//- k) + (s)\\ L 2 m ds. 
Jo 

Now, we use the Gronwall-Bellman lemma as in [3], Lemma A. 4, p. 156], and find that 

1/2 



4(v-k) + \\ 2 c0{[0tT] . H) + [ f \v(p-ky 

Jo Jn 





k 


< 










k 


< 









T 



\\Xk(t)\\ L 7/2 {n) \\d t p(t)\\ L 14/3 {n) dt 
—F II^PlL 7 /3(o,r;L 14 /3(n)) !|^fe||L7/*(0,T;iV2(m) • 

Next, we observe that the interpolation inequality (I2.20p (with p = 2, q = 6, r = 14/3, 
and = 1/7), together with the Sobolev inequality (12. 19ft . gives that 

IMUv3( 0> T;L«/3(n)) < {J ll v (*)llia^n) ll v (*)lli«(n) d *) 7 

-II II 1 / 7 f /" T |l mil 2 A 3/7 _ || II 1/7 || 1 1 6/7 

~~ W V Wl°°(0,T;H) y / H^WHi^)^ — IP Hl°°(0,T;J?) ll V lli2(0,T;L 6 (n)) 
*/ 

< II II 1 / 7 II II 6 / 7 

C \\ V \\l°°(0,T;H) W V Wl 2 (0,T;V)' 

and we denote by D the rightmost side of this inequality chain, when evaluated for 
v = dtp. We also remark that 



T 



1/2 -v 4/7 



X k\\L?/H0,T;LV*(Q)) = [j [J \ X k(x, t) \ 7/2 dxj dt 

.1/2 i.| g ^ 
Xk(x,t)\ dxj dt ? = ||'^fc||^2(o j T;L 4 (f2)) 



1/2 

\X k {x,t)\*dx] 

'o x Jn 

Hence, our estimate for (/i — k) + yields the following basic inequality: 

1 8/7 



fc)+f < ^i||^fc|| L 2 (0jT;L 4 (n)) for every k > p*, (4.6) 
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where D x = D / min{e, 1}, and where the norm |||-||| is defined by 

sup \\v(t)\\ 2 H + [ \Vv\ 2 for v E C°([0, T];H)n D 2 (0, T; V). 
te\o,T] Jo 



II 2 ._ 

te[o,T] 



We notice that the Sobolev inequality (I2.19P implies that 

\\v\\LZ(o,T;H(n)) < D 2 \\v\l for every v G C°([0,T]; H) fl D 2 (0,T; V), (4.7) 

where D 2 depends on Q and T, only. At this point, we select a sequence {kj} depending 
on a real parameter m > 1 as follows: 

kj := M(2 - 2 _i ) for j = 0,1,..., with M := mp*; (4.8) 

note that ko = M > Pq. Then, owing to (14.61) and (14. 7p . it is not difficult to check that 

- fy) ||Xfc i+1 11^(0,^(0)) < ~ %) + IU2(0,T;L4(n)) < ^K/i ~ kj) + \\ 

< kj D 1 D4X k f I !J i0jT ^ m . (4.9) 
Therefore, if we set 

Sj ■= ||^fc 3 HL 2 (0,T;L 4 (f7)) for j — 0, 1, ... , 

then the following inequality holds: 

< r ^ r D x D 2 Sf < W X D 2 VSf for j = 0, 1, . . . . 

Using [121 Lemma 5.6, p. 95], we conclude that Sj — > as j — > oo, provided that 

5 = 11X^11^(0,^(0)) < (4D lJ D 2 )- 7 2- 49 . (4.10) 

On the other hand, we notice that Xk = Xm, and we recall that M > p* and m = M//1q, 
by (14.8p . Moreover, we observe that Xm = 1 < (p — pl)/(M — p^) when /i > M, and that 
Xm = otherwise. Therefore, using (14.71) and (14. 6 j) with k = ko = M, we have: 

S ° - M - /ig ~ ^) + IU 2 (0,T;L4(C)) < j^Z^* 1^ ~ ^o) + i 

Dl Do ,, ii 8/7 -Di -Do , „, 1 8 1 8 



m 



1 ll / ^MLa(0 ) T;L4(n)) - _ 1 



We are now in a position to choose m :— 1 + D 1 D 2 |f2| 2//7 T 4 / 7 (4D 1 D 2 ) 7 2 49 . Then, m > 1 
and (I4.10p is satisfied. Consequently, 

||X2M||L 2 (o,T ; L 4 (n)) = hm Sj = 0, 

J-S-OO 

due to Beppo Levi's Monotone Convergence Theorem. This implies that p < 2M a.e. 
in Q , and the boundedness of p claimed in (12. 15[) is established. 

We are left with the task of proving that the limitations for p in (I2.15P do hold. We 
find it convenient to set: p* := inf^ p (recall that we assumed p* to be strictly positive). 
Moreover, we rewrite (I2.13P in the form: 

8d t p-Ap + f[(p)=g } where g:=p-f 2 {p), (4.11) 
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and we notice that g G L°°(Q), in view of the above proof and (12.31) . Consequently, in 
view also of (12 .4p . we can choose r* G (0, p*) such that /((r#) < g a.e. in Q. Then, we 
test (14.1 ip by — (p — r*) - and deduce that 

£ / |(p-r,)^(t)| 2 + /7 |V(p-r,)~| 2 - /7 (/i(p)-/J(r.))(p-r.)- 
z jo Jn Jo Jn 

= f I {fi(n)-g)( P -n)-<o. 

Jo Jn 

We conclude that (p — r*) - = and p > r* a.e. in Q. In a similar way, for a suitable 
r* < 1, we show that p < r* a.e. in Q by testing (14. lip by (p — 1 + r*) + . We conclude 
that solutions satisfy all of the requirements stated in (12.151) . □ 

5 Long-time behavior 

In this section, we prove Theorem 12.51 To this end, we fix any solution (p, p) to problem 
(I2.12p - (l2.14p . Our proof of the properties of the w-limit w(p, p) relies on a number of 
a priori estimates for (p, p), and on a well-known tool. For (p w ,p w ) any element of 
u)([i, p), and {t n } a corresponding time sequence of type (12.171) . we set 

H n (t) :=fi(t n + t), p n (t) := p(t n + t) fort>0, (5.1) 

and we study the sequence {(pn,p n )} on a fixed finite time interval [0, T]. Clearly, the 
pair (p n ,p n ) enjoys the same regularity as (p,p), and solves the equations 

(e + 2p n )d t p n + p n d t p n - Ap„ = a.e. in Q (5.2) 
Sd t p n - Ap n + f'(p n ) = p n a.e. in Q; (5.3) 

moreover, it satisfies the homogeneous Neumann boundary conditions and the Cauchy 
conditions: 

p n (0) = p(t n ) and p n (0) = p(t n ) a.e. in Q. (5.4) 
Our argument also relies on two basic identities, to be proved in the next lemma. 

Lemma 5.1. The following identities hold: 

8{d t pf - d t pAp + f'(p) d t p = ed t p + 2d t (pp) - Ap a.e. in Q; (5.5) 

$ [ t \\dtp(s)\\ 2 H ds + h\Vp(t)\\l+ [ f(p(t)) 
Jo A Jn 

= \ II Vpolli + / f(po) + e [ p(t) - e [ po + 2 / {pp){t) - 2 / p p , (5.6) 
z Jn Jn Jn Jn Jn 

for every t G [0, T] . 

Proof. We have from fl2~T2l that 

p d t p = Ap - ed t p - 2p d t p = Ap - ed t p - 2d t {pp) + 2d t p p . 
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By a simple rearrangement, we deduce that 

fi d t p = edtfi + 2d t (pp) - Ap. 

On the other hand, multiplication of ( 12 . 1 3[) by dtp yields: 

S(d t p) 2 - dtpAp + f'(p) dtp = pd tPl 

so that (15. 5p immediately follows by comparison. Next, identity (15. 6 p is arrived at by 
integrating (15 .5p over Q x (0, t) and by noting that, in view of the homogeneous Neumann 
boundary condition, Ap does not contribute to the integral. □ 

We proceed with proving some a priori estimates. In so doing, we depart from our 
general rule, and write c for constants that do not depend on the final time T, although 
they are allowed to depend on the element (/i^,^) of the w-limit under consideration. 
Whenever a dependence of c on the parameter T cannot be excluded, we stress this 
possibility by writing ct- Morever, without any loss of generality, we assume that e < 1. 

First a priori estimate. Just as we did for (I3.18p . we immediately deduce that 

f\\Vp(s)\\ 2 H ds+ S -\\p(t)\\ 2 H + [( P fi 2 )(t)<c foreveryt>0. (5.7) 
Jo z Jn 

This implies, in particular, that 

r+oo 

/ \\V p(t)\\ 2 H dt < +oo and ||/!„||l-(o,t ; h) < c. (5.8) 
Jo 

Second a priori estimate. We recall (15.61) and estimate some terms of its right-hand 
side. We have: 



e / »{t) + 2 / (pp)(t) Ke 1 ' 2 / p{t) + 2 / {p l l 2 p){t) <2\V\+ e\\p{t)\\ 2 H + / {pp 2 ){t). 
Jn Jn Jn Jn Jn 

On the other hand, (15.71) holds and / is bounded from below. Hence, (15. 6p yields: 

r* i 

S \\dtp(s)\\ 2 H ds+-\\Vp(t)\\ 2 H <c for every t > 0, (5.9) 
Jo 2 



whence we have that 



r+oo 

/ \\9tP(t)\\ 2 H dt < +oo and \\p n \\L°°(o,T;V) < c- (5.10) 
Jo 

Third a priori estimate. We formally test (12.131) by — Ap and integrate over Q x 
[t n ,t n + t). Due to the convexity of fi and the boundedness of f% (compare with the 
derivation of (I3.16P ). we get 

z Jt n Jn 

X rtn+t f 1 ftn+t f 1 rtn+t f 

< ^p(tn)\\ 2 H + C Wp\ 2 + -J \^P\ 2 + -J / H 2 - 

J tji J Q J tn J Yl J tn J f2 
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We note that the first term on the right-hand side is bounded, since p(t n ) is weakly 
convergent to p w in V . Hence, owing to (15. 9 \ and (15 .7p . we can conclude that 

rtn+t r 

l|Vp(tn +t)f H + / / |Ap| 2 < c T for every t G [0,T] and every n. 

Jt n JSl 

By comparison with ( 12. 131) . and by exploiting elliptic regularity, we deduce that 

rtn+t rt-n+t 

/ \\f[(p)\\ 2 H dt+ \\p(t)\\ 2 wdt<c T . 

In terms of p n , all this reads: 

\\pn\\L2(o,T;W) < c T and || f[ (p n ) \\ L *(o,T;H) < or- (5.11) 



Fourth a priori estimate. By rewriting (15. 2p in the form 

p n 8 t p n Ap n 



dtPr, 



e + 2p n e + 2p n ' 



(5.12) 



and owing to the homogeneous Neumann boundary condition, we are entitled to write 
the following equation in V* in the framework of the Hilbert triplet (V, H, V*): 



d t p{t) v 



Pn{t) d t pn(t) V 

e + 2p n (t) 



[ Vp n (t) ■ V- 
Jn t 



2Pn{t) ' 



(5.13) 



for a.a. t G (0,T) and for every v G V. Starting from this equation, we can prove 
a bound for d t p n in L p (0,T; V*) for some p > 1. For a while, we argue for a.a. t G 
(0,T) and estimate each term on the right-hand side separately. Owing to the Sobolev 
inequality (12.191) . we get 



p n {t) d t p n (t) v 



< C\\p n 



< e 1 \\Vn{t)\\LHn) \\dtp(t)\\ L 2 {n) \\v\\ L 4 {Q) 



2p n (t) 

L*{n) \\dtp(t)\\ L 2(n) \\v\\v 



(5.14) 



On the other hand, using the interpolation inequality f l2.20j) and the Sobolev inequality 
once more, we obtain that 



^ \\ AAII 3 / 4 II z'+MI 1 / 4 II C+MI 3 / 4 II ^Ml 1 / 4 

L\n) < \\Pn{t)\\ L6{n) \\Pn{t)\\ L2{u) < C\\p n {t)\\ v ||/i n ' + 



II 



Consequently, by accounting for the L°° bound of (15. 8p . we derive from (15. 14j) that 

Pn(t) d t p n {t) V 



n e + 2p n {t) 



<c\\Pnit)\\T\\dtp{t)\\H\\v\\ V . 



(5.15) 



This takes care of the first addendum in on the right-hand side of (15.131) . As to the second, 
we have: 

' < \\p n (t)\\v \\\7((e + 2 Pn (t)y 1 v)\\ H , (5.16) 



/ V/i„(t)-V 



e + 2p ri {t) 
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where it remains for us to estimate the last norm. By applying the Leibniz rule and 
making use of the Holder interpolation and the Sobolev inequalities, we get: 

\\V((e + 2 Pn (t))~ 1 v)\\ H < 2e- 2 \\vV Pn (t)\\ H + e- 1 \\Vv\\ H 

< c(||v||l*(o) ||Vp B (t)|| L 4 (n) + \\Vv\\ H ) < c(\\V Pn (t)\\ L 4 {n) + l)||v|| v 

< c(\\V Pn (t)\\%\ n) ||Vp n (t)||^ + < C(\\p n (t)\\% 4 \\ Pn (t)\\\! 4 + l)\\v\\ V . 

Hence, on accounting for the L°° bound of (I5.10p . we see that (I5.16P becomes: 

J Vp n (t) ■ V g - ^ | < c\\p n (t)\\v(\\Pn(t)\\% 4 + 1) Hv ■ (5.17) 
Since v G V is arbitrary, by combining (15.131) . (I5.15p . and (15.171) . we arrive at: 

\\d t fin(t)\\v* < C{\\p n mf \\d tP {t)\\ H + U n {t)\\ V \\p n (t)fw + (5.18) 

for a. a. t G (0, T) and for every n; by estimating each term on the right-hand side, we are 
going to find some p > 1 such that 

\\d t ^n\\LP(0,T;V*) < C T - (5.19) 

Now, due to (I5.8p . the last term is bounded in L 2 (0,T). As to the first, we observe that 
the functions 

t^\W(t)\\T and t ^ \\d t p(t)\\ H 

are bounded, respectively, in L 8//3 (0, T) by ( 15. 8ft . and in L 2 (0,T) by (15. 10ft ; hence, their 
product is bounded in L 8//7 (0,T), by the Holder inequality. Finally, the middle term 
of (15. 181) can be treated in a similar way, with the use of (15.81) and the first inequality 
in (15.111) . Therefore, (I5.19P does hold, with p = 8/7. 

Conclusion. Our estimates (15.81) . (15.101) . and (I5.19P ensure that (fi, p) is a bounded 
and weakly continuous (H, y)-valued function. Hence, the first part of the statement in 
Theorem 12.51 follows from the general theory (see, e.g., [HI p. 12]). We pass to the study 
of the cu-limit. 

Recalling (15. lip and using standard weak and weak star compactness results, we see 
that there is a triplet (p,oo, Poo, (fioo) such that 

Pn ->> Poo weakly star in L°°(0, T; H) n L 2 (0, T; V), (5.20) 
Pn-^Poo weakly star in H^O.T; H)r)L°°(0,T;V) D L 2 (0,T;W), (5.21) 
f'APn) -)• 0oo weakly in L 2 (0, T; H), (5.22) 

at least for some subsequence. Our first aim is to prove that /Xqo is a nonnegative constant, 
i.e., that ^(Xjt) = p, s for a.a. (x,t) G Q for some p, s G [0, +oo); and, to prove that p^ 
is time independent, i.e., that poo(i) = p s for a.a. t G (0,T) for some p s G W. Secondly, 
we want to prove that the pair (fi s ,p s ) found in such a way is indeed a steady state and 
coincides with the given pair (/j, 

U! 1 Pui ) ■ 

From the first bounds of (15. 8p and (15.101) . we immediately deduce that 
|V/i n | -> and d t p n -> strongly in L 2 (0, T; H). 
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This implies that /i^ is space independent and is time independent. Thus, we can 
write Poo(t) = p s for a. a. t E (0, T), for some p s E W. Moreover, (I5.2ip implies strong 
convergence: 

Pn^ Poo strongly in C°([0, T]; iJ)flL 2 (0, T; V) (5.23) 

(see, e.g., [TTJ Sect. 8, Cor. 4]). Therefore, f^Pn) converges to /^(poo), e.g., strongly in 
L 2 (0,T; H), and it is clear that 

- Ap^ + ^ = poo - /a(Poo) a.e. in Q. (5.24) 

Collecting f!5.23j) and (I5.22p . and recalling Lemma 1.3, p. 42, in p], we conclude that 
< Poo < 1 and 0oo = /{(poo) a.e. in Q. Therefore, (15.241) becomes: 

< p s < 1 and - Ap s + f{p s ) = Poo a.e. in Q, 

and we deduce that p^ is time independent as well. Thus, p 00 (a ; , £) = Ps for a. a. (x, t) E Q 
for some constant p s . Furthermore, p s is nonnegative, since p n > for every n. This 
concludes the proof that (p s ,p s ) is a steady state. 

It remains for us to show that (p s ,p s ) coincides with (p^p^). From (I5.23P we see 
that p n (0) converges strongly in to poo(0) = p s ; on the other hand, p n (0) = p(t n ) 
converges weakly in V to p w , by assumption; hence, p s = p w . A similar argument holds 
for p s and p w , because p n converges strongly in C°([0, T]; V*). Indeed, p n , is bounded 
in L°°(0,T;H), by ([53]); on the other hand, fl5TT9|) holds with p = 8/7 > 1, and the 
embedding if C V^* is compact; hence, the desired convergence follows from [TTl Sect. 8, 
Cor. 4]. This completes the proof of Theorem 12.51 
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